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Abstract 

Aims. We model the present-day number and properties of ultracompact X-ray binaries (UCXBs) in the Galactic Bulge. The main 
objective is to compare the results to the known UCXB population as well as to data from the Galactic Bulge Survey, in order to learn 
about the formation of UCXBs and their evolution, such as the onset of mass transfer and late-time behavior. 

Methods. The binary population synthesis code SeBa and detailed stellar evolutionary tracks have been used to model the UCXB 
population in the Bulge. The luminosity behavior of UCXBs has been predicted using long-term X-ray observations of the known 
UCXBs as well as the thermal-viscous disk instability model. 

Results. In our model, the majority of UCXBs initially have a helium burning star donor. Of the white dwarf donors, most have helium 
composition. In the absence of a mechanism that destroys old UCXBs, we predict (0.2 - 1.9) x 10 5 UCXBs in the Galactic Bulge, 
depending on assumptions, mostly at orbital periods longer than 60 min (a large number of long-period systems also follows from the 
observed short-period UCXB population). About 5-50 UCXBs should be brighter than 10 35 erg s -1 , mostly persistent sources with 
orbital periods shorter than about 30 min and with degenerate helium and carbon-oxygen donors. This is about one order of magnitude 
more than the observed number of (probably) three. 

Conclusions. This overprediction of short-period UCXBs by roughly one order of magnitude implies that fewer systems are formed, 
or that a super-Eddington mass transfer rate is more difficult to survive than we assumed. The very small number of observed long- 
period UCXBs with respect to short-period UCXBs, the surprisingly high luminosity of the observed UCXBs with orbital periods 
around 50 min, and the properties of the PSR J1719-1438 system all point to much faster UCXB evolution than expected from angular 
momentum loss via gravitational wave radiation alone. Old UCXBs, if they still exist, probably have orbital periods longer than 2 h 
and have become very faint due to either reduced accretion or quiescence, or have become detached. UCXBs are promising candidate 
progenitors of isolated millisecond radio pulsars. 
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1. Introduction 

Ultracompact X-ray binaries (UCXBs) are low-mass X-ray bina- 
ries with observed orbital periods shorter than ~ 1 h, indicating 
a compact, hydrogen deficient donor star ( |Vila|197T| Paczyriski 
1 1 98 1 1 [Sien kiewicz 1984). The donor overflows its Roche lobe, 
and lost matter is partially accreted by a neutron star or black 
hole companion. Because of the compact orbit, mass transfer 
is driven by orbital angular momentum loss via gravitational 
wave radiation (e.g. Kraft et al.|1962 Paczyriski| 1 967| |Faulkner 



l97Tj |Pringle & Webbink||1975| |Tutukov & Yungelson||1979) . 
For an overview of the relevance of studying UCXBs, see e.g. 
|Nelemans & Jonker| ( |20T0l i. 

The Galactic Bulge is a suitable environment to dis- 
cover UCXBs because of the high local star concentration. 
Furthermore, the Bulge is an old stellar population that contains 
few young X-ray sources such as high-mass X-ray binaries and 
core-collapse supernova remnants. Due to their rarity, UCXBs 
are not found nearby in the Galaxy. Also, while observable in 
X-rays, they are too faint to be identified in other galaxies (e.g. 
|Voss & Gilfanov|2007D . 



The Galactic Bulge Survey (GBS) ( |Jonker et al.|201 1| > is an 
X-ray and optical survey focused on two 6° X 1 ° regions centered 
1 .5° to the North and South of the Galactic Center. One of the 
goals of the GBS is to investigate the properties of populations 
of X-ray binaries in order to constrain their formation scenarios, 
especially the common-envelope phase(s). 

The present study aims to predict and explain GBS results 
regarding the number and luminosities of UCXBs by means 
of binary population synthesis and UCXB evolutionary tracks, 
thereby contributing to a better understanding of the formation 
and evolution of UCXBs. We also compare our results to the or- 
bital periods and chemical compositions of the known UCXB 
population. In Sect. [2 we describe our assumptions on the star 
formation, stellar and binary evolution, and the observable char- 
acteristics of evolved UCXBs. The results follow in Sect. [3] 
where we present the modeled present-day population and its 
observational properties. In Sect. |4l we compare our resu lts with 
the population synthesis studies by |Belczynski & Taam (2004a), 
|Zhu et aL| l |2012a) , and |Zhu et al.|p012b| ), as well as to obser- 
vations, and we discuss various implications. We conclude in 
SectE 
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2. Method 

The study of the evolution of the UCXB population consists of 
several steps. First, the star formation history of the Galactic 
Bulge and the binary initial mass function prescribe which types 
of zero-age main sequence binaries form in the Bulge, and 
when they form. The binary population synthesis code SeBa 



et al. 2012) 1 is used to simulate the evolution of this popu 
lation of zero-age main sequence binaries. During the evolu- 
tion of the population, all UCXB progenitors are selected at a 
certain moment after the supernova explosion that leaves be- 
hind the neutron star or black hole. More detailed evolutionary 
tracks are used to trace the subsequent evolution. This yields 
the present-day number and intrinsic parameters of the UCXBs 
in the Galactic Bulge. Finally, using long-term observations by 
the Rossi X-ray Timing Explorer All-Sky Monitor (RXTE ASM) 
( |Bradt et aL]|1993t |Levine et aI71|1996| > as well as the accre- 
tion disk instability model, the modeled UCXB parameters are 
translated into observational parameters using the results of van 
Haaften et al. ( 2012c| l. The result is a prediction of the present- 
day observable population. 

2.1. Population synthesis 

The binary population synthesis code SeBa models the evolu- 
tionary transformations of a population of binary stars based on 
a distribution of initial binary parameters. It follows the evo- 



lution of stellar components using analytic formulas by Hurley 
et al. ( 2000[ l, taking into account circularization due to tidal in 



teraction, magnetic braking, gravitational wave radiation, mass 
exchange via Roche-lobe overflow, common envelopes, and em- 
pirical parameterizations of wind mass loss. For a more exten- 
sive description of SeBa we refer to Portegies Zwart & Verbunt 



( fl^9^ , |Neiemans et al.lpOOT] ), and |Toonen et al.| ( |2012| i 

As is common, we use parameterizations to describe the 
common-envelope process. In this study, the density profile of 
the donor envelope is parametrized by A = 1 /2, the efficiency 
with which orbital energy is used on unbinding the common en- 
velope by a?cE - 4 (justified by explosive shell burning in mas- 
sive stars during the common-envelope phase, Podsiadlowski 
et al. 2010 1 and the specific angular momentum of the enve- 
lope after it has left the system, relative to the specific angular 
momentum of the pre-common-envelope binary, by y = 7/4. 
In choosing a value of ctce^ = 2 for massive stars we follow 



Portegies Zwart & Yungelson (1998); Yungelson et al. (2006); 
Yungelson & Lasotafi2008|l. The y- and acE-prescriptions are 



used as describ ed in 
has been used (Zoccali et al. 12003 



bonen et al.|([20T2] >. A metallicity Z = 0.02 

"IT" 



The metallicity is held con- 
stant in the population synthesis simulations as well as the subse- 
quent tracks because of the relatively short episode of star forma- 
tion (Sect. |2~2] ). For the kicks acquired by nascent neutron stars 
we use the velocity distribution suggested by Paczyhski ( 1990) 
with a dispersion of 270 km s . Because these parameters are 
very uncertain, in Sect. |4] we will consider the effect of varying 
the common-envelope parameters and the supernova kick veloc- 
ity distribution from our standard values. 



2.2. Star formation history and initial binary parameters 

The star formation history of the Galactic Bulge can be approx 
imated by a Gaussian distribution with a mean fi — -10 Gyr and 
a standard deviation <x = 0.5 - 2.5 G yr, where the total mass 
of stars that are formed is 1 x 10 10 M Q ( |Clarkson & Rich|2009 ' 
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( Portegies Zwa rt & Verbunt|1996||"Nelemans et al.|2001| |Toonen| ~ 
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Figure 1. Star formation history of the Galactic Bulge as a 
Gaussian distribution for mean fi = -10 Gyr and two values 
of the standard deviation cr. The total mass of the stars formed is 
1 x 10 10 Mq. Time = corresponds to the present. Star formation 
is assumed to start at Time = -13 Gyr. 



Wyse 2009). Star formation is assumed to start 13 Gyr before 
present. For a narrow distribution the star formation is concen- 
trated around 10 Gyr in the past, but the cr — 2.5 Gyr distribution 
has an important tail of recent star formation. In this paper we 
consider cr = 0.5 Gyr and cr = 2.5 Gyr (Fig. [TJ, representing 
a relatively instantaneous burst of star formation, and additional 
recent star formation, respectively. 

We find the initial binary parameters by primary-constrained 
pairing ( TKouwenhoven et "al7||2008| >. We derive the initial mass 
function for the primary components by combining the initial 
mass function for stellar systems by Kroupa (2001 ) with an esti- 
mate of the binary fraction as a function of mass. The mass ratio 
of the secondary and the primary is drawn from a flat distribu- 
tion between and 1. The 1 x 10 10 M Q of stellar mass is found to 
contain 5.2 x 10 9 binary systems. We simulated 1 million bina- 
ries with a lower primary mass limit of 4 M Q (because systems 
with a lower primary mass do not produce a supernova event in 
either component), and another 4.3 million binaries with a lower 
primary mass limit of 8 M Q , after it became clear that systems 
with lower primary masses do not produce UCXBs. Using the 
binary initial mass function we calculated to how many binaries 
in the full mass range of primaries (0.08 - 100 M Q ) this simula- 
tion corresponds. The resulting population was then multiplied 
by a factor (of 14.45) to scale to the entire Bulge population. 

For an analytic derivation of the binary initial mass function 
we refer to Appendix [A] 

2.3. Formation scenarios 

We consider three UCXB-progenitor classes, each defined by the 
stellar type of the donor at the time it fills, or will fill, its Roche 
lobe: 

Class 1. White dwarf with a neutron star or black hole compan- 



Iben et al. 1995 Yungelson 



ion (Tutukov & Yungelson 1993 
let al.|2002D , 

Class 2. Helium burning star with a neutron star or black 
hole companion ( jSavonije et al. 1986| Iben & Tutukov 1987 
Yungelson|20"08) , 



Class 3. Evolved main sequence star of about 1 M© with a neu- 
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tron star or black hole companion (fTutukov et al.|1985[|Ne lson 



et al.|1986||Fedorova & Ergma|19 89; Pyly ser & Savonije| 1989 



Podsiadlowski et al.||2002| |Nelson &R appaport 2003 1 |van der 
Sluys et al.|2005a||Lin et al.|201 1) . 



These classes include all the binary systems that may even- 
tually evolve into an UCXB (Belczynski & Taam 2004a van der 



|Sluys et al.|20 05a; Neleman s et al.|2010) > - in some models in 
volving accretion-induced collapse of a white dwarf or a neutron 
star, the donor star has already transferred mass before the for- 
mation of the eventual accretor, a neutron star or a black hole, 
respectively. The detailed tracks follow the evolution of the he- 
lium burning donor starting immediately after its formation, and 
the main sequence donor immediately after the supernova event. 
The white dwarf donor tracks start at the onset of Roche-lobe 
overflow. In each of the detailed tracks the mass transfer is con- 
servative as long as the mass transfer rate does not exceed the 
Eddington limit. If the mass transfer is faster than that, accre- 
tion at the Eddington limit is assumed, and the mass that is lost 
from the system carries the specific angular momentum of the 
accretor. 

The initial system parameters (component masses and or- 
bital periods) and major events during the evolution towards an 
UCXB are described below for each class. 

2.3.1 . White dwarf donor systems 

The evolution of UCXBs with a white dwarf donor can be di- 
vided into two main categories based on whether the primary 
(initially more massive star) or secondary component becomes 
a supernova. A supernova explosion of the secondary star is 
possible when it gains mass by hydrogen accretion from the 
primary (e.g. Tutukov & Yungelson 1993; Portegies Zwart & 



Verbunt|1996J|van Kerkwijk & Kulkarni|1999[|Portegies Zwart 
& Yungelson||1999||Tauris & Sennels||2000| >. A supernova ex- 
plosion of the secondary probably never produces a black hole, 
neither does the primary turn into a helium white dwarf after the 
supernova, because it starts out too massive. Thus, all secondary- 
supernova systems have a carbon-oxygen or oxygen-neon white 
dwarf donor and a neutron star accretor. Because the high stel- 
lar mass required for a supernova explosion of the primary is 
relatively rare due to the steep initial mass function, in our simu- 
lations a significant fraction of the systems (13% of the carbon- 
oxygen white dwarf systems and 36% of the oxygen-neon white 
dwarf systems) experience their supernova in the secondary star. 
Systems with a black hole accretor are rare, about 0.2% of all 
white dwarf systems. All black holes form from the primary and 
have a ~0.6 - 0.8 M Q carbon-oxygen white dwarf companion. 

Carbon-oxygen white dwarf system^] are 1 .5 times more 
prevalent than oxygen-neon white dwarf systems in our simu- 
lations, and ~ 30 times more prevalent than helium white dwarf 
systems (combining primary and secondary supernovae). 

Supernova explosion of the primary This category can be sub- 
divided by the predominant white dwarf composition: helium, 
carbon-oxygen or oxygen-neon. 

Evolution starts with a zero-age main sequence binary in 
which the primary is a massive star (M > 8 M if the sec- 
ondary is to become a helium or carbon-oxygen white dwarf, 
and M > 10 M Q in the case of an oxygen-neon white dwarf 
companion) that evolves off the main sequence first. Systems 



1 The envelopes of some white dwarfs contain helium, but we refer 
to them as carbon-oxygen white dwarfs because the envelope will be 
lost relatively early. 



that eventually produce a helium white dwarf donor have ini- 
tial orbital periods ranging mainly from 1 to 100 yr. For sys- 
tems that produce a carbon-oxygen or oxygen-neon white dwarf 
donor the orbital periods lie mostly between 0.1 and 1000 yr. In 
the case of neutron star accretors, the primary expands during 
the Hertzsprung gap or as a giant and fills its Roche lobe, fol- 
lowed by mass transfer to the companion. The secondary cannot 
accrete all of this mass and is engulfed in a common envelope 
( |PaczynskT 1976 1. The envelope is expelled before the two com- 
ponents merge and the exposed helium core and the main se- 
quence secondary are left behind with an orbital separation sev- 
eral tens of times smaller than before. After a few 10 Myr, the 
helium star turns into a giant (which may lead to a subsequent 
phase of Roche-lobe overflow) and explodes as a core-collapse 
supernova, leaving behind a neutron star. Black hole progenitors 
are Wolf-Rayet stars, which lose a large fraction of their mass 
before evolving into a helium giant, and then collapse to form a 
black hole. 

Secondary stars that eventually become a helium white dwarf 
donor had a zero-age main sequence mass between 1.4 and 
2.3 M , whereas for carbon-oxygen white dwarf donors this 
range is 2.3-7 M , where 2.3 M G is the maximum mass of single 
stars that undergo the helium flash. A small fraction started with 
a higher initial mass. Progenitors of oxygen-neon white dwarf 
secondaries have a mass between 7 and 1 1 M on the zero-age 
main sequence and do not become a supernova due to severe 
mass loss (e.g. Gil-Pons & Garcia-Berro|2001 1. (There is some 
overlap with the progenitor mass range of carbon-oxygen white 
dwarfs - the end product depends on whether burning stops be- 
fore of after carbon ignition.) After the supernova explosion has 
occurred, the secondary evolves off the main sequence. As a 
subgiant, it initiates a common envelope with the neutron star, 
shrinking the orbit by another factor of a few tens. The core 
cools into a helium white dwarf (< 0.35 M Q ) or, after a helium 
burning and helium giant stage, a carbon-oxygen white dwarf 
(0.35 - 1.1 M ), or even, after carbon burning, an oxygen-neon 
white dwarf (> 1.1 M , |Gil-Pons & Garria-Berro|2001| >. Orbital 
angular momentum loss via gravitational wave radiation further 
shrinks the orbit until the white dwarf eventually overfills its 
Roche lobe, which happens at an orbital period of a few min- 
utes. 



Supernova explosion of the secondary In this scenario, the to- 
tal binary mass needs to be at least 9 M if the primary becomes 
a carbon-oxygen white dwarf and 12 M if the primary becomes 
an oxygen-neon white dwarf. The primary transfers several so- 
lar masses to the secondary in a stable manner (avoiding a com- 
mon envelope - initially the secondary must have a mass of at 
least 0.55 times the primary mass) while ascending the red giant 
branch (e.g. |Tauris &~S ennels 2000). Eventually the core be- 
comes a helium star, or a carbon-oxygen or oxygen-neon white 
dwarf. The secondary, which becomes the more massive compo- 
nent of the system, evolves off the main sequence and initiates 
a common envelope. The orbit shrinks, and 30 - 70 Myr after 
the binary formation the secondary explodes as a supernova and 
produces a neutron star. In systems which remain bound after the 
supernova explosion, the primary will eventually reach Roche- 
lobe overflow as a relatively massive (> 0.7 M ) carbon-oxygen 
white dwarf or an oxygen-neon white dwarf (> 1.1 M ). The 
relatively high initial mass of the primary precludes less mas- 
sive carbon-oxygen white dwarfs. The initial orbital period in 
this scenario can be much shorter than in the primary-supernova 
scenarios, down to a few days. The initial stellar masses lie be- 
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tween 4.5 and ~ 10 M for the primary and 4-9 M for the 
secondary (if the former becomes a carbon-oxygen white dwarf) 
and about 5-12 M for both the primary and secondary (if the 
former becomes an oxygen-neon white dwarf). 

The onset of mass transfer from the white dwarf Most white 
dwarf-neutron star systems merge upon the onset of mass trans- 
fer. For a 1 .4 M neutron star companion, white dwarfs with a 
mass higher than ~ 0.83 M G experience dynamically unstable 
mass transfer, assuming a zero-temperature (completely degen- 
erate) mass-radius relation for the donor (e.g. Yungelson et al 
2002 van Haaften et al.||2012b[ ). This assumption implies that 



these white dwarfs have cooled considerably by the time they 
eventually fill their Roche lobe, although tidal heating and ir- 
radiation before the onset of mass transfer may counteract this 
for a short timej^This leads to runaway mass loss on the dy- 
namical timescale of the donor, followed by accretion of part 
of the disrupted white dwarf via a disk around the neutron star 
(see e.g. van den Heuvel & Bonserria||1984| [Fryer et al.| |l999 



Paschalidi s et al.|201 1| Metzger|2012 1. Furthermore, in systems 
with a donor mass larger than ~0.38 M Q ( |Yungelson et al.|2 002 



van Haaften et al.|2012b[l (this value is only weakly sensitive to 



accretor mass) the accretor will be unable to eject enough trans- 
ferred matter from the system by isotropic re-emission, where 
most arriving matter leaves the vicinity of the accretor in a fast, 
isotropic wind powered by accretion ([Soberman et al. 1997 



Tauris & Savonije 1999). This also leads to a merger. Therefore, 



systems that are unstable due to either a dynamical instability 
or insufficient isotropic re-emission have been removed from the 
sample. These instabilities only occur in systems with a white 
dwarf donor, because of the negative donor mass-radius expo- 
nent and the small donor size (hence, small orbit) at the onset 
of mass transfer. Dynamical instabilities may occur in systems 
with helium or main sequence donors if they have masses con- 



siderably higher than considered in this study (see e.g. Pols & 
|Marinus|1994[ >. 

In our simulation 97.4% of all white dwarf systems have a 
donor with a mass higher than 0.38 M and do not survive the 
onset of mass transfer. This includes 99.1% of carbon-oxygen 
(solid line in Fig. |2]> and all oxygen-neon (dashed line in Fig. |2j 
white dwarf systems. In about 80% of the surviving white dwarf 
donor systems, the donor is a helium white dwarf (dotted line in 
Fig. [2}, in the remainder it is a carbon-oxygen white dwarf. All 
surviving systems experienced the supernova explosion in the 
primary star and host a neutron star. 

If a white dwarf donor with a mass higher than the 0.38 M 
isotropic re-emission limit has a non-degenerate surface layer, 
the system may not merge immediately upon the onset of mass 
transfer, but it will once this layer has been lost. 

Two thirds of the white dwarf donor systems start to trans- 
fer mass to the neutron star within 2 Gyr, but some systems take 
much longer (Fig. [3j. This is the case for all white dwarf types. 
White dwarfs can take very long to start mass transfer depending 
on the width of the initial orbit, since the orbital decay of bina- 
ries consisting of a neutron star and a white dwarf is caused by 
gravitational wave radiation only. 

Evolutionary tracks For each donor composition, the evolution 
after the stage described in Sect. |2.3.1| follows the tracks de- 



2 If the white dwarfs are warm, their radius is larger at a given mass 
(Bilds ten|2002)|Deloye & Bildsten|2003| > and therefore the orbital peri- 
ods we find are strictly lower limits. 
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White dwarf mass at onset Roche-lobe overflow (M s ) 

Figure 2. Total time-integrated number of UCXBs with a helium 
(dotted line), carbon-oxygen (solid line) or oxygen-neon (dashed 
line) white dwarf donor at the onset of mass transfer to a neutron 
star (including merging systems). 
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Figure 3. Delay time distribution between the zero-age main 
sequence (ZAMS) and the onset of mass transfer to a neutron 
star for UCXBs with a white dwarf (excluding merging systems, 
solid line), helium star (dashed line) and main sequence (dotted 
line) donor. 



scribed in van Haa ften et al.| ( |2012b| l. Initially, the white dwarf 
donor has not yet cooled and therefore is larger than a zero- 
temperature white dwarf of the same mass. While the donor 
loses mass, its radius is held constant until it equals the zero- 
temperature radius of the same mass (this is justified by the rapid 
mass loss the donor initially experiences). From this point on, the 
zero-temperature radius (Zapolsky & Salpeter 1969; Rappaport 
|et al.|1987] > is used, which increases with further mass loss. The 
initial neutron star mass is taken to be 1.4 M and its radius 12 
km HGuillot et aLl[20TT] |Steiner et afJ|2l)T3] l. The evolution of 
UCXBs with degenerate donor stars is governed by angular mo- 
mentum loss through gravitational wave radiation, which forces 
mass transfer via Roche-lobe overflow. 
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2.3.3. Main sequence donor systems 

Main sequence donors have mostly evolved from 1.0- 1.2 M Q 
secondaries that started transferring mass after orbital decay due 
to magnetic braking ( |van der Sluys et a l.|2005a I. After the mag- 
netic field disappears (because the star becomes fully convective 
as a result of mass loss), gravitational wave radiation becomes 
the dominating angular momentum loss mechanism, continuing 
the orbital shrinking. In this scenario, the initial period and donor 
mass need to fall within relatively narrow ranges in order to suf- 
ficiently evolve the main sequence star. Moreover, the magnetic 
braking must be sufficien tly efficient ([va n der Sluys et al. 2005b ) 
which it probably is not (Qu eloz et al.|1998] >. Depending on the 
extent of hydrogen depletion in the stellar center, systems can 
reach a minimum orbital period between 10 and 80 min, where 
~ 80 min is the lower limit for hydrogen-rich donors (Paczyhski 
[T9811 



Figure 4. Sample of helium star-neutron star UCXB tracks, for 
initial helium star mass range 0.35 -1.0 M Q and initial orbital 
period range 20-200 min by Nelemans et al. (2010 1 (solid lines), 



and our white dwarf donor track extensions (dashed lines). 



We have used stellar evolutionary tracks by van der Sluys 
|et al.| ( |2005a| l for binaries with an 0.7 - 1 .5 M Q main sequence 
star and a 1 .4 M neutron star at initial orbital periods between 
0.50 and 2.75 days. The donor metallicity Z — 0.01. These tracks 
describe the orbital period and mass transfer rate as a function 
of time, as well as the core and surface compositions. 



2.3.2. Helium burning donor systems 



2.4. Behavior of old UCXBs 



The supernova in low-mass helium burning star systems occurs 
in the primary, which has an initial mass M > 8 M Q . Most sys- 
tems start out with an orbital period between 0.1 and 100 yr. 
Most of the helium stars form less than 500 Myr after the zero- 
age main sequence when the secondaries, with initial masses of 
2.3 - 5 M Q experience hydrogen shell burning (Savonije et al 



1986) and lose their hydrogen envelope in case B mass trans- 
fer ( Kippenhah n & Weigert|[T967 1. They fill their Roche lobes 



within another ~ 200 Myr, which is much earlier than UCXBs 
with white dwarf donors. In part this is due to the requirement 
that the helium star has not yet turned into a white dwarf before 
the onset of mass transfer (thereby disqualifying itself from the 
helium burning donor sample) which constrains the size of the 
initial orbit. We do not find systems with a black hole accretor. 

We have used stellar evolutionary tracks for systems with 
an 0.35 - 1 .0 M Q helium star and a 1 .4 M neu tron star at ini- 
tial orbital periods^] between 20 and 200 min (Nelemans et al. 
|2010) table in electronic article), part of which is shown in Fig. [4] 
These tracks were made in the same way as the tracks for sys- 
tems with white dwarf accretors in Yungelson (2008 ). The donor 
metallicity Z = 0.02. Because the donors at the end of the 
tracks are degenerate, we have extended these tracks by making 
a smooth transition to the zero-temperature white dwarf evolu- 
tion. The tracks describe the orbital period, mass transfer rate, 
donor mass, and core and surface compositions as a function of 
time. In Fig.|4] the helium stars live up to 400 Myr. After the on- 
set of mass transfer (vertical part of the tracks), the orbits shrink 
until the period minimums, then expand towards the bottom right 
of the figure. For each individual 'zero-age' UCXB system pro- 
duced by SeBa, the track that best matches its donor mass and 
orbital period has been used. 



3 The initial orbital period is the period directly after completion 
of the common-envelope phase that leaves behind the helium star 
( Yungelson|2008) . 



Figure |4] suggests that once the donor has become degenerate, 
UCXBs 'uneventfully' reach long orbital periods and very low 
mass transfer rates. This is probably not the case. Instead, at low 
mass transfer rate a thermal-viscous instability in the accretion 
disk ( |Osaki|[T9"74] |Lasota|[200T] l can cause UCXBs with a suf- 



ficiently low mass transfer rate to become transient (Deloye & 
Bildsten 2003 ). This usually implies that these systems are vis- 
ible only during outbursts when the disk is in a hot state, which 
is only a small fraction of the time, and not during the quies- 
cent state when the disk is cold and gaining mass. Furthermore, 
due to accretion of angular momentum, a neutron star accretor 
in an U CXB can be recycled to a spin period between one and a 
few ms (|Bisnov atyi-Ko gan & Komberg|1974 Alpar et aL| l982 



Radhakrishnan & Srinivasan 1982). Combined with a low mass 
transfer rate, the magnetosphere may transfer angular momen- 
tum from the neutron star to the accretion disk, thereby accel- 
erating orbiting disk matter and counteracting accretion, known 
as the 'propeller effect' ( Davidson & Ostriker||1973| |Illarionov| 
|& Sunyaev||1975[ ). As a result, the inner accretion disk, source 
of most X-ray radiation, can become disrupted by the magneto- 



sphere. See van Haaften et al. (2012b) for more details on the 
thermal-viscous disk instability and propeller effect in UCXBs. 
Finally, at low donor mass, high-energy radiation from the neu- 
tron star, the magnetosphere and the accretion disk may evapo- 



rate the donor, or detach it from its Roche lobe ( Kluzniak et al 



|1988[ |van den Heuvel & van Paradljs 1988 Ruderman et al 



1989a; Rasio et al. 1989). Hot, low-mass donors may suffer from 
a dynamical instability caused by a minimum value of their mass 
in the case of a constant core temperature (Bild sten|2002[ ). 

Each of these mechanisms can potentially diminish the visi- 
bility of UCXBs. Because it is impossible to precisely quantify 
at which stage of the evolution (if at all) these mechanisms be- 
come important, and to what degree, we do not remove UCXBs 
from the sample, instead we will discuss the implications on the 
population of old UCXBs in the Discussion (Sect. |4j). 
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2.5. Present-day population 

The present-day number and system parameters of UCXBs in 
the Galactic Bulge can be found by evaluating the evolutionary 
stage of all simulated systems at the present time. The most in- 
teresting parameters are the orbital period, mass transfer rate and 
surface composition, because these can be inferred from obser- 
vations. The orbital periods of observed systems can be found 
from periodic modulations in the light curve or spectrum, al- 
though this is usually very difficult for UCXBs (e.g. Nelemans 
|et al.|2006| |in't Zand et al.||2007| ). Since the transferred matter 
originates from the surface of the donor, the occurrence and rel- 
ative abundance of elements in the donor can be inferred from 
X-ray ( |Schulz et al.|2001|>, ul traviolet ( |Homer et al.|2002) , and 
optical (Nelem ans et al.|2004)> spectr a, and more indirectly, type- 
I X-ray bursts ( in't Zand et al. 2005). The mass transfer rate can- 
not be directly determined observationally. However, because 
the energy output of an X-ray binary is for a large part provided 
by the gravitational energy release of the accreted matter, the 
mass transfer rate strongly influences the luminosity of the sys- 
tem, which can be observed. 



2.5.1. Bolometric luminosity 
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Figure 5. UCXB variability: fraction of time that a source emits 
above a given luminosity for 14 UCXBs, including two candi- 
dates with tentative orbital periods, adapted from van Haaften 
|et al.| ( |2012c| ). The numbers associated with the curves indicate 
the orbital periods in minutes. 



We employ two methods of converting the modeled mass trans- 
fer rate to bolometric luminosity. 

An observational method is to match the modeled systems to 
real systems and assume that the modeled system behaves simi- 
lar to the real system in terms of emission. We match a modeled 
UCXB to the real UCXB with the nearest orbital period. The 
relevant parameter of the emission behavior is the fraction of 
the time a source radiates at a given bolometric luminosity, mea- 
sured over a sufficiently long timespan. We use 16-year obser- 
vations by the RXTE ASM to determine this behavior for the 14 
known UCXBs (including two candidates) for which ASM data 
is available. Figure[5]shows this behavior for sources when they 
are visible well above the noise level ( |van Ha aften et al. 2 012c| l. 
ASM X-ray luminosity was converted to bolometric luminosity 



using an estimate by |in't Zand et alT (2007 1. At a given time, 
the luminosity of an UCXB is randomly drawn from either the 
individual ASM observations that make up this time-luminosity 
curve, or (most of the time) from the faint-end extrapolations of 
the curves in Fig.|5l( |van Haaften et al.|2012c| l. These extrapola- 
tions are constructed in such a way that the average luminosity 
of the luminosity distribution is equal to the time-averaged lu- 
minosity of the source as observed by the ASM. The amount of 
time that a given source spends at a particular luminosity trans- 
lates into the number of sources in a population at the same lu- 
minosity. 

The second method, of a more theoretical nature, is to con- 
vert a system's modeled mass transfer rate to luminosity, using 
predictions by the disk instability model (Sect. |2.4| i in the case 
of long-period UCXBs. According to this model the mass trans- 
fer rate must exceed a critical value in order to be stable and 
the source to be persistent, i.e., visible at a relatively high lu- 
minosity (almost) all the time. A crude estimate for the critical 



mass transfer rate in the case of an irradiated disk is given by in't 
Zandet al.|(|2007), b ased on |Dubus et aI] ( [l999) ; |Lasota|p001) ; 
Menouet al.| l |2002) 



5.3 x 10" 11 / 
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with M a the accretor mass, P or b the orbital period, and / is a 
factor accounting for the disk composition; f ~ I for carbon- 
oxygen disks and f ~ 6 for helium disks. 

When the time-averaged mass transfer rate exceeds the criti- 
cal value, the bolometric luminosity L is assumed to be constant 
at 

_ GMaMa 

with G the gravitational constant and R d the accretor radius. 

Sources with a time-averaged mass transfer rate below the 
critical value are assumed to be visible only during outburst 
stages. The predictions by the thermal-viscous disk instability 
model regarding the degree of variability of sources is supported 
by observations ( |van Paradijs|1996"{ |Ramsay et al.|2012) |Coriat 
|et al.|2012) . The duty cycle (fraction of the time the source is in 
outburst) is 

DC = -^2-, (3) 

^outburst 

where L avg is the time-averaged luminosity based on the theo- 
retical mass transfer rate, (Eq. |2|, and ^outburst is the luminosity 
during outburst, derived by Lasot a~et al.| ( |2008) 



^outburst ~ 3.5 X 10 



37 



/P \1- 67 

I h ) 



erg s 



(4) 



which is consistent with observations of outbursts in UCXBs 
(e.g. |Wu et aL||2010) . The period of this cycle is not relevant 
here. We neglect the decay in the light curve after an outburst. 
Furthermore, we do not predict the luminosity of systems that 
are in quiescence, which in fact has been assumed to be zero in 
the above method. 

Both methods have advantages and shortcomings. The ASM 
observations have a rather high lower limit in converted bolo- 
metric luminosity, of ~ 10 37 erg s" 1 at 8.3 kpc, the estimated 
distance to the Galactic Center (Gille ssen et al.|[2009) >. Below 
this luminosity, we have to rely on extrapolations. Also, variabil- 
ity on a timescale longer than 16 yr cannot have been observed 
by the ASM. On the other hand, the data do not strongly rely 
on uncertainties in models, as is the case for the disk instability 
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model method. The ASM data show that UCXBs with a simi- 
lar orbital period can behave rather differently, for instance XTE 
J1751-305 (42.4 min orbital period) and XTE J0929-314 (43.6 
min), or 4U 0513^0 (17.0 min) and 2S 0918-549 (17.4 min). 
Equation ([TJ, however, shows that in the disk instability model 
the critical mass transfer rate and hence luminosity is largely de- 
termined by the orbital period and composition. A longer orbital 
period does not automatically imply a lower luminosity, as evi- 
denced by M 15 X-2 (20.6 min orbital period) and 4U 1916-05 
(49.5 min). 

In general, the ASM data show that a clear distinction be- 
tween persistent and transient behavior is not justified (van 
Haaf ten et al.|2012c) . Almost all systems are visible above the 



ASM detection limit only sporadically. Still, short orbital period 
systems are typically visible more often at a given luminosity, 
i.e., they have a (slightly) higher time-averaged luminosity. Even 
though the available sample is small, individual unusual behav- 
ior is expected to partially cancel out for the population as a 
whole, because some modeled UCXBs will be matched to a real 
UCXB that is brighter than typical for its orbital period, while 
others will be matched to one that is fainter than typical. 

2.5.2. Donor composition 

The donor surface compositions of the modeled UCXBs are pre- 
dicted using the helium-star donor and main sequence donor 
tracks, as well as the white dwarf types from the population syn- 
thesis model. These predicted compositions can be compared 
to observations of real systems, in the Bulge and elsewhere. 
Donors that start mass transfer as a white dwarf can be helium 
and carbon-oxygen white dwarfs (Sect. 2.3.1 1. In the latter case 



we assume 30% carbon and 70% oxygen by mass, based on 
the most common eventual compositions in the helium burning 
donors (Sect. |3.4) , Donors that start mass transfer as a helium 
burning star can also produce helium-carbon-oxygen donors due 
to an interrupted helium burning stage. For the subsequent tracks 
for these systems, we use the mass-radius relation for degen- 
erate donors composed of a mixture of 60% helium, 30% car- 
bon and 10% oxygen, a choice based on the dominant tracks by 
|Nelemans et al . (2010) as will be discussed in Sect. |3.4| We note 
that the degenerate tracks are not very sensitive to the composi- 
tion (as long as there is no hydrogen), so these simplifications are 
justified. Matter processed in the CNO cycle has a high nitrogen- 
to-carbon abundance ratio, whereas helium burning converts this 
to a low ratio. Consequently, the nitrogen-to-oxygen ratio is a 
good test for the formation channel because it can discriminate 
between a history as a helium white dwarf donor or a helium 
burning donor ( Nelemans et al. 2010|l. 



Based on the overview in 



van Haaften et al.| ( |2012c| >, the 



compositions of observed UCXBs can be summarized as be 
ing roughly equally distributed over helium and carbon-oxygen 
compositions. There is no clear dependency on the orbital pe- 
riod, although helium composition may be more common among 
systems with a long orbital period (> 40 min) ( |van Haaften et al. 
2012c). The surface composition of very low-mass donors cor- 



responds to the (inner) core composition of the object before it 
started transferring mass. 

3. Results 

3.1. Birth rates and total number of systems 

Convolving the star formation history (Fig. [TJ with the delay 
times of the onset of mass transfer (Fig. [3} yields the birth rate 




-12 -10 -8 -6 
Time (Gyr) 

Figure 6. Birth rate of systems reaching Roche-lobe overflow 
against time for UCXBs from the white dwarf donor channel 
(solid lines), helium-star donor channel (dashed lines) and main 
sequence donor channel (dotted lines). Time = corresponds to 
the present. Black lines correspond to a star formation history 
width <x = 0.5 Gyr, gray lines to cr = 2.5 Gyr. 



distributions, shown in Fig. |6]for burst-like (black) and extended 
(gray) star formation epochs. Except for part of the white dwarf 
channel and the main sequence channel, the delay times are 



much shorter than the age of the Bulge (Sect. 2.2 1. In the case of 
cr = 0.5 Gyr, 98% of the white dwarf donor systems have started 
Roche -lobe overflow before the present, whereas 100% of the 
helium-star donor and 84% of the main sequence donor systems 
have. In the case of cr = 2.5, these percentages are 97%, 100% 
and 77%, respectively. The main distinguishing feature between 
the three classes (Sect. [23) is the most recent time at which mass 
transfer can begin. Initially wide systems from the white dwarf 
donor channel can still start mass transfer at the present, whereas 
main sequence donor systems and especially helium-star donor 
systems cannot, unless they have formed relatively recently (star 
formation history width cr = 2.5 Gyr). The rate of helium burn- 
ing donor systems reaching Roche-lobe overflow closely follows 
the star formation history. 

Upon the onset of mass loss, the donor radius increases im- 
mediately for fully degenerate white dwarf donors, and after ap- 
proximately 100 Myr for helium burning donors, once the donor 
has become sufficiently degenerate following the extinction of 
nuclear fusion (this happens some time after the period mini- 
mum). The orbital period decreases in the case of helium burn- 
ing or main sequence donors, whereas the period increases with 
mass loss for systems with degenerate donors (in each channel). 
If main sequence donor systems become ultracompact, this typ- 
ically happens ~ 3 Gyr after the onset of mass transfer, and mass 
transfer starts after 2-6 Gyr after the formation of the binary 
( |van der Sluys et al.|2003a) . 

The total number of UCXBs with a white dwarf or helium- 
star donor, shown as the solid and dashed lines in Fig. [7] initially 
follows the star formation rate and later on approaches an upper 
limit as star formation slows down. The number of systems be- 
low a given orbital period initially resembles the instantaneous 
star formation rate more closely for short orbital periods. After 
the peak in star formation rate, the number of systems below 
a given period keeps increasing as long as more new systems 
form than old systems are removed from the given sample due 
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Figure 7. Number of UCXBs from the white dwarf donor chan- 
nel (solid lines), helium-star donor channel (dashed lines) and 
main sequence donor channel (dotted lines). Time = corre- 
sponds to the present. For each line style, the three lines rep- 
resent the full population (upper) and the systems with orbital 
periods shorter than 60 min (middle) and shorter than 20 min 
(lower). A star formation history width cr = 0.5 Gyr has been 
used. 
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Figure 8. Present-day orbital period distribution for UCXBs 
from the white dwarf donor channel (solid lines), helium-star 
donor channel (dashed lines) and main sequence donor channel 
(dotted lines). The elements next to the lines indicate the most 
abundant element(s) at the surface of the donor, hence in the 
transferred matter. The star formation history has width cr = 0.5 
Gyr. 



to their increasing orbital periods. The numbers of systems from 
the three different classes decline at different rates corresponding 
to their respective recent birth rates (Fig. [6]!. 

3.2. Present-day population 

While the evolution of the population is interesting in itself, the 
population today can be used to validate the results. In the case 
of a star formation history distribution width <x = 0.5 Gyr, shown 
in Fig. [8] most systems at the present are old, and have expanded 
to an orbital period of ~ 80 min. Since evolution slows down at 
longer periods, systems tend to 'pile up' j^] Differences in donor 
composition lead to different present-day orbital periods. This is 
the case even among hydrogen-deficient compositions because 
during most of the evolution, the donor mass is low enough for 
Coulomb physics to be important to the stellar structure, or even 
dominating degeneracy pressure. Coulomb interactions cause a 
donor that is composed of 'heavy' elements such as carbon and 
oxygen to have a smaller radius than donors with lighter compo- 
sition, such as helium, of the same mass (Zapo lsky & Salpeter| 
1 1969) 1. A larger donor radius (at each mass) results in a longer 
orbital period at each mass, but also at each age (because less 
time is spent at a given orbital period)]^] 

For cr = 0.5 Gyr, all UCXBs with orbital periods shorter 
than 1 h started Roche-lobe overflow from a white dwarf donor 
(Fig. [8]), long after the formation of the system. Most of these 
systems host a helium white dwarf. The main sequence channel 
contributes a negligible number of UCXBs and can only be dis- 
tinguishable (in principle) via donor compositions. For cr = 2.5 
Gyr, shown in Fig. [9] there is also recent star formation. This 
produces a population of young UCXBs that descended from 
helium burning donor systems (or still have a helium burning 
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Figure 9. Same as Fig. [s] except cr = 2.5 Gyr. 



donor), with orbital periods shorter than 1 h. The steep cut-off 
at the long-period end of several curves is due to the assumption 
that star formation suddenly starts 1 3 Gyr before present. 

Combining all donor compositions, the result is a current 
~ 1.9 x 10 5 population of UCXBs, mostly at long orbital pe- 
riod (60 - 90 min). The total number of UCXBs in each class is 
3.5 x 10 4 (18%) with white dwarf donors, 1.56 x 10 5 (81%) with 
helium-star donors, and 5.1 x 10 2 (0.3%) with main sequence 
donors. The number of modeled systems with orbital periods 
shorter than 60 min is 1.5 x 10 3 for cr = 0.5 Gyr, and 7.4 x 10 3 
for cr = 2.5 Gyr (0.8% and 3.8% of the population, respectively). 
We note that these numbers are rather sensitive to assumptions 
in the model, and could be lower by an order of magnitude, as 
will be discussed in Sect. [4] 
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Figure 10. Present-day luminosity distribution of UCXBs in 
the Bulge based on Rossi XTE All-Sky Monitor observations, 
after incorporating the accelerated evolution of the systems 
(Sect. 3.3.1 1. The star formation history has width cr = 0.5 Gyr 



(solid line) and cr = 2.5 Gyr (dashed line). 



3.3. Observable population 



As described in Sect. |2.5.1[ in order to determine what we can 
observe at high luminosity we have to convert the modeled pop- 
ulation to luminosities. 



3.3.1. RXTE All-Sky Monitor 

In the first method, we apply the observations of known UCXBs 
by the RXTE ASM (Fig. |5]_Sect. lOTT) to the modeled popula- 
tion (Figs. [8] and [9] Sect. 133. Modeled UCXBs with an orbital 



period longer than 60 min are left out because of the absence of 
known real systems with such periods (i.e., they are assumed 
never to reach luminosities above ~ 10 34 erg s~'). The time- 
averaged luminosity of most UCXBs with orbital periods longer 
than 40 min is approximately two orders of magnitude higher 
than expected from the gravitational-wave model ( van Haaften 
et al.|2012c i. This implies that either the observed sources are 
atypically bright, or that they show normal behavior, but evolve 
much faster than if driven only by gravitational wave losses (the 



implications will be discussed in Sect. 4.2 1. In each case, due 



to energy conservation, we need to reduce the number of bright 
sources at each orbital period by a factor that corresponds to the 
ratio between the gravitational-wave luminosity and the actual 
observed luminosity, given by van Haaften et al. (2012c their 
Fig. 3). Figure 10 shows the resulting number of bright UCXBs 



predicted by the ASM data. 

For star formation history width cr = 0.5 Gyr, about 40 sys- 
tems are expected to be visible as bright sources at a given time. 
For cr = 2.5 Gyr, this number is larger, ~80, because recent star 
formation causes more young systems to exist, which have not 
yet reached orbital periods of 60 min. The cut-off at the faint 
end of the histogram is an artifact of the assumed linear extrap- 
olation to the faint behavior. This also results in a relatively high 
minimum luminosity. In reality, especially sources with orbital 
periods longer than ~ 40 min are expected to be very faint (i.e., 
much fainter than suggested by a linear extrapolation) at least 
some fraction of the time, which means the cumulative luminos- 
ity distribution flattens at faint luminosities, causing a tail at the 



Figure 11. Present-day luminosity distribution of UCXBs in the 
Bulge based on the disk instability model (Sect. 2.5.1 1. Line 



color distinguishes between helium (thick black lines), carbon- 
oxygen (thick gray and thin black lines) and helium-carbon- 
oxygen (thin gray lines) donor compositions. The thick lines cor- 
respond to the white dwarf donor channel, the thin lines to the 
helium burning channel. The star formation history has width 
cr = 0.5 Gyr (solid lines) and cr = 2.5 Gyr (dashed lines). 



low-luminosity end of Fig. 10 The distribution decreases with 



increasing luminosity in a similar way as the luminosity distri- 
bution of the representative individual observed UCXBs (Fig. [5}. 

3.3.2. Disk instability model 

The second method relies on converting theoretical mass 
transfer rates to luminosities using the disk instability model 
(Sect. 2.5.1 1, the result of which is shown in Figs. 1 1 and 12 The 



total number of bright (> 10 35 erg s~') sources (either persistent 
or in outburst) is 34 for cr = 0.5 Gyr and 51 for cr = 2.5 Gyr. 
The vast majority of these are persistent (short-period) sources, 
and therefore the number is larger in the case of recent star for- 
mation. For the same reason, the white dwarf donor channel 
(Sect. ~~ 



2.3.1 



dominates this population, and most should have 
helium or carbon-oxygen donors. The helium burning channel 
is expected to contribute at most a couple of bright sources, in 
outburst at long orbital period (60 - 80 min). 
Figure 



11 



shows the luminosity distribution of bright 
sources. UCXBs with a luminosity between 
are the ones with the shortest orbital periods, below 



10 36-37 erg s -l 



20 min. 

From here, the number of sources at a given luminosity in- 
creases towards fainter luminosities because these sources have 
longer period derivatives and lower time-averaged mass transfer 
rates, until a sharp cut-off defined by the longest orbital period 
at which sources are still considered persistent. Carbon-oxygen 
dominated sources are persistent to longer periods and lower lu- 
minosities because accretion disks composed of carbon-oxygen 
are more stable than helium dominated disks, see Eq. ([T}. For 
cr = 0.5 Gyr the number of short-period, persistent, systems is 
negligible. At longer periods (> 40 min), the duty cycle (Eq.[3]l 
determines the number of sources in outburst. The duty cycle 
decreases below 10 at orbital periods longer than 60 - 70 min, 
depending on donor composition. During their rare outbursts, 
sources are temporarily bright at ~ 10 37 ~ 38 ergs -1 , as follows 
from Eq. Q. Their peak luminosities are used in Fig. [TT| 
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Figure 12. Orbital period distribution of the predicted bright 
population of UCXBs in the Bulge at the present based on the 



disk instability model. For details see Fig. 1 1 



Two peaks can be distinguished in the lines representing 
the white dwarf channel in Fig. 11 (thick lines). The peak at 



~ 10 erg s consists of systems with orbital periods just longer 
than the critical period because these still have a relatively high 
duty cycle. The second peak at ~ 10 37 6 erg s _1 consists of long- 
period systems because these are very numerous and distributed 
over a relatively narrow interval of orbital periods. The duty cy- 
cle at a given orbital period is higher for helium burning systems 
owing to their larger size and, because their average density is 
set by the orbital period, correspondingly larger mass. Hence, 
their time-averaged mass transfer rate at a given orbital period is 
also higher. 

In Fig. 12 the orbital period distribution is shown for the 
same population of bright sources as in Fig. [TT] The jumps of 
these distributions correspond to the respective low-luminosity 



ends of the distribution in Fig. 1 1 The cut-off period of persis- 
tent sources (at 30-40 min) lies at a longer orbital period for the 
systems with a helium burning donor origin compared with sys- 
tems with a white dwarf origin. The reason is that these donors 
have a higher temperature than originally white dwarf donors, 
and therefore the time-averaged mass transfer rate is higher at 
the same period. This causes the disk to remain stable (and the 
sources to be persistent) up to a longer period. Again we see 
that carbon-oxygen donor systems are persistent up to a longer 
orbital period than helium-dominated donor systems. Transient 
systems with orbital periods longer than ~ 40 min are rarely in 
outburst and at most a handful have a high luminosity at a given 
time. 



3.4. Donor surface composition 

The helium-star donors have partially turned into carbon-oxygen 
white dwarfs during their evolution, depending on their mass 
and evolutionary stage at the onset of Roche-lobe overflow (de- 
termined by the initial mass and orbital period). When the star 
starts mass transfer after filling its Roche lobe, burning is ex- 
tinguished quickly (Sav onije et al.|[l986| l, and at this stage the 
core mass fraction of helium varies between a few percent to al- 
most 100% ( |Nelemans et al.|2 010). Figure 13 shows the surface 
abundances at the present day, assuming a narrow star forma- 
tion history (<x = 0.5 Gyr). Two thirds of the systems end up 
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Figure 13. Surface abundances of helium versus carbon (black 
squares) and oxygen (white triangles) for donors in the helium- 
star channel at the present time in the case of a sufficiently nar- 
row star formation history. These correspond to the core abun- 
dances at the end of the tracks by Nelemans et al. (2010). The 
surface area of a symbol is proportional to the number of systems 
in the corresponding track. 



with less than 10% helium on their surface. Systems that started 
out with a short orbital period generally have a higher helium 
mass fraction, because these had less time to burn helium be- 
fore the onset of mass transfer. The abundances depend on the 
temperature at which helium and carbon burning takes place. A 
higher temperature causes a higher helium burning rate, produc- 
ing more carbon. Later, the carbon abundance reduces in favor 
of oxygen. The scatter in Fig. 13 is therefore due to differences 



in core burning temperature caused by different stellar masses. 

UCXBs produced via the white dwarf donor channel have 
donors mostly composed of either helium material (produced 
in the CNO cycle), or carbon-oxygen. The ratio between both 
types is about 3:1 but strongly depends on the efficiency of 
isotropic re-emission, which strongly affects the number of 
carbon-oxygen white dwarf donor systems that survive the onset 
of mass transfer (Sect. |2~3~T| i. 

Evolved main sequence donors with an initial mass of 1.0 - 
1.2 M Q reach a helium surface abundance Y * 0.9 once they 
become ultracompact, the remaining part being hydrogen and a 
metallicity Z = 0.01, which has been present from the start ( |van| 
der Sluys et al.|2005a| l. 

Because the core material is exposed early on (e.g., at an 
orbital period shorter than ~ 20 min for the helium burning sys- 
tems, Nelema ns et al.|2010 1, and the core is homogeneous due 
to convection during its burning stages, the chemical composi- 
tion is expected not to change with increasing orbital period, and 
therefore the same for the total and the observable population. 

3.5. Collective emission as function of orbital period 

Even though variability behavior determines the number of 
UCXBs in outburst and their luminosity, the collective luminos- 
ity of all UCXBs at a given orbital period is in principle not de- 
pendent on variability, because the time-averaged mass transfer 
rate of an UCXB is a relatively straightforward function of or- 
bital period. During the evolution of an UCXB, its evolutionary 
timescale increases with age and orbital period. This means that 
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Figure 14. Energy released per unit interval of orbital period for 
an UCXB with a zero-temperature helium white dwarf donor 
and a 1.4 M Q neutron star accretor or a 10 M Q black hole ac- 
cretor. Note that the quantity on the vertical axis should not be 
interpreted as a luminosity; the time unit represents change in 
orbital period rather than passing time. 

there exist many more systems at longer orbital period. On the 
other hand the time-averaged mass transfer rate decreases with 
age and orbital period, so a long-period source has a lower time- 
averaged luminosity. The total energy output of a source, or of 
the population as a whole (under the assumption of a constant 
star formation rate), at a given orbital period is an indication of 
at which orbital periods systems are likely to be observed. 

The amount of energy emitted by an UCXB per unit orbital 
period is given by 



dE 



dp 



orb 



L 

Kb 



GM a dM d 
~^dP^ 



(5) 



where E is the emitted energy and the donor mass. This re- 
lation is illustrated in Fig. [14] The donor mass decreases much 
faster at shorter orbital periods (Ma 



Pj-^ (van Haaften et al. 



2012b|), i.e., dM d /dP orb cc p^ 3 ), and since the donor mass is the 



luel lor the luminosity, systems emit much more energy at short 
orbital periods, not only per time interval (their luminosity) but 
also per orbital period interval. For instance, an UCXB will emit 
~ 12 times as much energy during its evolution from 20 to 21 min 
as it does between 60 and 61 min. The consequence is that the 
short-period systems dominate the collective X-ray output of an 
UCXB population, unless the star formation rate decreases very 
fast. Depending on the variability of systems and the sensitivity 
of the instrument used, this could very well result in short-period 
systems dominating the visible population^] 

4. Discussion 

We predict ~ 1 .9 x 10 s of UCXBs in the Galactic Bulge, predom- 
inantly at orbital periods of > 70 min, but also a few thousand 
systems with orbital periods shorter than 60 min (but mostly 
longer than 40 min). Based on RXTE ASM observations, about 



6 Not to be confused with the known population, which is composed 
by individual observations at different times. In the known population, 
long-period systems will eventually dominate provided they can be seen 
during outbursts 



Table 1. Size of the modeled UCXB population in the Galactic 
Bulge for different model parameters. 





Kick distribution 


Number of UCXBs (x 10 4 ) 


2 


Paczyhski 


19 


2 


Maxwellian 


14 


0.2 


Paczyhski 


4 


0.2 


Maxwellian 


2 



40 - 80 of these sub-hour UCXBs should be visible at high lu- 



minosities of > 10 erg s (Fig. 10 1, depending on the star 
formation history. Also, ~ 35 - 50 bright UCXBs with orbital 
periods < 30 min (i.e., persistent) should be visible above such 
a luminosity based on grav itational energy release and the disk 
instability model (Figs. |l l|and 12 1 ^\ 

The combined common-envelope parameter oqe^ for mas- 
sive stars may be lower than the value of 2 that we used ( Voss 
& Tauris 2003). Decreasing this value to 0.2, as well as us- 
ing a Maxwellian kick velocity distribution with a dispersion 
of 450 km s _1 , rather than the distribution by |Paczynski ( 1990} , 
would reduce the number of UCXBs formed by a factor of ~ 8, 
as fewer systems will survive the common-envelope stage or the 
supernova explosion. In that case, the total number of UCXBs 
in the Bulge we predict is ~ 2 x 10 4 , and the number of bright 
systems becomes ~5 - 10. Table[T]shows the number of UCXBs 
in our model for various combinations of common-envelope ef- 
ficiency and neutron star kick velocity distribution. Furthermore, 
the slope of the initial mass function at high stellar mass is also 
uncertain. A steeper slope (resulting from earlier studies such as 
|Kroupa et al.| [T993 ) leads to a smaller fraction of massive stars 
and therefore fewer UCXBs. A different choice for the initial 
component mass pairing may also reduce the number of UCXBs 
by an order of magnitude (Belczynski & Taam 2004a). 

We can distinguish several disagreements with observa- 
tions. First, no UCXBs with orbital periods longer than 60 min 
have been discovered, faint or bright, in the Bulge or else- 
where. Second, no bright UCXBs with a short orbital period 
(< 30 min) have been identified in the Galactic Bulge. Third, 
only three UCXBs with orbital periods between 40 and 55 
min, XTE J1807-294 (|Markwardt et al.|2003] ), XT EJ1751-305 
(Mar kwardtetal.|2002} and SWIFT J 1756.9-2508 ( |Krimm et al.| 
|2007[ ), are presumably located in the Bulge, based on their posi- 
tions in the sky, as their distances are not known. 

As for the predicted ~ 1.9 x 10 s long-period systems, the 
probable existence of three observed UCXBs with orbital peri- 
ods shorter than 55 min in the Bulge can be used to calibrate the 
formation rate of UCXBs, independent of population synthesis. 
This yields a much larger number of UCXBs than three for sys- 
tems with a period longer than 55 min, based only on the rapid 
increase of the evolutionary timescale (set by gravitational wave 
radiation) with orbital period. For instance, UCXBs are expected 
to reach an orbital period of 55 min within ~ 2 Gyr after the onset 
of mass transfer ( |van Haaften et al.|2012b| ). 

Even though the disk instability model predicts these long- 
period systems to be in outburst so rarely that few are expected 
to be bright, in quiescence they would still be detectable by sen- 



7 The observed UCXBs with orbital periods around 20 min are not 
clearly detected by ASM most of the time, even though they are pre- 
dicted to be persistent by the disk instability model. Incorporating the 
fact that they are bright less than 100% of the time, the predicted number 
of bright systems becomes lower, depending on the precise luminosity 
distribution of these systems. 
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sitive instruments at ~ 10 31 33 erg s 1 (e.g. Bildsten & Rutledge 



200Tj|Heinke et al.|2003||Belczynski & faam|2004b| ). 



The three observed UCXBs that are located in the direction 
of the Bulge have undetermined donor compositions, though 
SWIFT J1756.9-2508 is thought to have helium composition 
( |Krimm et al.| [2007). This cannot be used to constrain the star 
formation history, as UCXBs with helium white dwarf donors 
can form with a wide range of delay times. 



in the Disk, Klypin et al.||2002 )plThe overprediction of bright, 
persistent UCXBs is therefore not unique to our study. 

Recently, Zhu et al. ( 2012b| > performed a population synthe- 
sis study of Galactic UCXBs with neutron star accretors, and 
predicted 5-10 x 10 3 systems in the Galaxy, depending on 
neutron star birth kicks. As in our study, the helium burning 
donor channel was the most common. Notable differences with 
our work are that these authors found a large number of UCXBs 
with a carbon-oxygen white dwarf origin, and a peak in the or- 
bital period distribution near 40 min. 



4.1. Comparison with previous studies 



4.2. Overprediction of UCXBs with long orbital period 



|Belczynski & Taam| ( |200 4a) performed a population synthesis 
study of (primordial) UCXBs in the Galactic disk. The main dif- 
ference between their and our results is that they found a total 
of 478 UCXBs with orbital periods shorter than 80 min in the 
disk at the present epoch, about three orders of magnitude fewer 
than our result, per unit star forming mass. It is not clear what 
causes this discrepancy, although these authors use different ini- 
tial binary parameters than we do (Appendix lAj, for example a 
steeper high-mass slope in the initial mass function that leads to 
fewer massive stars, relatively. Also, their initial primary masses 
leading to UCXBs span a narrower range. Of all UCXBs in their 
simulation, about 20% have a black hole accretor, all of which 
form via the accretion-induced collapse of a neutron star. This 
percentage strongly depends on the assumed upper mass of a 
neutron star (2 M ), the mass retention efficiency of the accret- 
ing neutron star and the evolutionary stage at which the common 
envelope happens. As in our study, these authors did not find any 
(surviving) UCXBs with a black hole that was formed directly in 
the collapse of a massive star. They found that about 90% of the 
neutron star accretors form via the accretion-induced collapse of 
an oxygen-neon-magnesium white dwarf, a scenario our model 
does not produce. In our simulations, 81% of the UCXBs start 
mass transfer from a helium burning donor, compared to 40% in 



Belczynski & Taam ( 2004a I for the Galactic field. This is not un- 



reasonable given the uncertainties in e.g. the onset of mass trans- 
fer from a white dwarf, and differences in assumptions between 
both studies. The number of persistent sources predicted by the 
disk instability model depends sensitively on the orbital period 
separating the persistent and transient sources, because most of 
the predicted persistent sources have orbital periods only slightly 
shorter than this critical period. Details in accretion disk models 
(e.g. X-ray irradiation) and composition can make a large dif- 
ference. We predict about 0.02% of the UCXBs to be persistent 
(Sect. 3.3.2|l, a much smaller fraction than found by Belczynski 



& Taam (2004a) (2.2%), but their number applies to the Galactic 



disk, which has a several orders of magnitude higher ongoing 
star formation rate, and therefore more young UCXBs, which 
have short orbital periods. 

Our number of persistently bright UCXBs (35 - 50 from the 
disk instability model, using our standard parameters) can also 
be compared with the number of persistent UCXBs with white 
dwarf donors (600 - 900) predicted by |Zhu et al!] ( |2012a| l for the 
whole Galaxy, if one takes into account the difference in adopted 
cut-off donor mass for persistent behavior. These authors found 
that UCXBs with donor masses lower than 0.03 M Q are transient 
whereas our limit lies around 0.02 M . Using their limit, our es- 
timate would reduce to roughly ten, which scales within a factor 
of a few with their number given the stellar mass ratio between 
Bulge (1 x 10 10 M ) and Galaxy (an additional 4 - 6 x 10 10 M Q 



An important clue towards what may happen at long periods 
comes from the long-term ASM data. The reason for the dif- 
ference in predictions by the ASM and disk instability model 
(Sect. [3T3J lies in the ASM observations that the UCXBs with 
orbital periods 40 - 55 min are approximately two orders of 
magnitude brighter than theoretically expected from the time- 
averaged mass transfer rate, assuming mass transfer is driven 
exclusively by gravitatio nal wave radiation in a binary with a 
(semi-)degenerate donor ( van Haaften et al.|2012 c )p] Assuming 
that the observed systems have been displaying normal behav- 
ior during the 16 years of RXTE observations, additional angular 
momentum loss besides that from gravitational wave radiation 
would cause a higher mass transfer rate at the same orbital pe- 
riod, and therefore a higher time-averaged luminosity (see also 
Ruderma n et al.|1989b[ ). As mentioned in Sect. |2.4| an efficient 
physical mechanism for additional loss of angular momentum 
from the system is a wind from the donor, induced by irradi- 
ation from the accretion disk or millisecond pulsar. In black 
widow systems, which host a millisecond pulsar and a low-mass 



(< 0.2 M Q ) companion in a < 10 h orbit ( |King et aL]|2005 
such donor evaporation has been observed (e.g. Fruchter et al 
1988). This scenario has also been proposed to be happening to 
the unusually light (~ 10~ 3 M G ) detached companion to the mil- 
lisecond pulsar PSR J 17 19-1438 ( |Bailes et al.pOTT]) via either 
the white dwarf or helium burning donor channels (|van Haaften 



et al. 2012a) or the evolved main sequence donor channel (with 
an orbital period minimum at ~ 45 min, Benvenuto et al.|2012| ), 
though the latter scenario does not produce a carbon-oxygen rich 
donor. 

The recently discovered spin-powered millisecond gamma- 
ray pulsar PSR J13 1 1-3430 system dPletsch et al.|2012b, with an 
orbital period of 93.8 min ( |Romani|2012[ |Kataoka et al.||2012| ) 



and an evaporating helium donor ( |Romani et al.|2012[ ) "supports 
the hypothesis that UCXB evolution is strongly influenced by 
donor evaporation. Given the low hydrogen abundance, donor 
evaporation, orbital period, pulsar spin period, and minimum 
companion mass, PSR J131 1-3430 could very well be an UCXB 
descendant on its way to becoming a millisecond radio pulsar 
system like PSR J1719-1438. 

The non-detection of UCXBs with periods longer than ~ 60 
min, when the donor is still expected to be much more massive 
than the companion in PSR J1719-1438, suggests the existence 
of another mechanism that hides UCXBs with low donor masses 
and low mass transfer rates (the 60 min limit is uncertain due 



to the small observed sample). The propeller effect (Sect. 2.4 1 



Zhu et al. 



(2012a) bund a much larger number of UCXBs with 



helium burning donor than our study, but this is expected given the large 
difference in recent star formation between Bulge and Disk. 

9 These ASM luminosities are direct observations, unrelated to the 
extrapolations of the curves shown in Fig.|3]and discussed in Sect. 3.3 
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is the most promising mechanism to explain this non-detection. 
The rotational energy of the millisecond pulsar is sufficient to 
make long-period UCXBs with very low mass transfer rates 
(~ 1(T 13 M Q yr -1 ) much fainter, since it causes arriving matter 
to be unbound ( van Haaften et aL|2012b) . The propeller effect 
could still allow for a (very) low rate of accretion that would 
prevent radio emission and make the sources visible in the ul- 
traviolet (owing to their low disk temperatures and possibly dis- 
turbed inner accretion disks). Furthermore, radio emission from 
a millisecond pulsar itself, once switched on after an interrup- 
tion in mass transfer, is capable of preventing accretion (Burderi 
etal.|2001l|Fu&Li|2011) 



Using the Chandra X-Ray Observatory (Weisskopf et al. 
2002), the Galactic Bulge Sur vey has found 1234 X-ray sources 
in 8.3 deg 2 ( Jonker et al. 201 1\ so far, most of which have not yet 
been identified. Although many are expected to be foreground 
Cataclysmic Variables or non-ultracompact X-ray binaries, this 
number of systems found in approximately 5% of the total area 
of the Bulge on the sky is at least consistent with a large pop- 
ulation of faint X-ray binaries. Also, a potentially large popula- 
tion of sub-luminous X-ray transients with neutron star accretors 
exists near the Galactic Center (|Sakano et al.||2005| |Wijn ands 
eraLl[2006l |Degenaar & Wijnands|2009| |2010> . These systems 
have (intrinsic) peak luminosities near ~ ltF" 5 erg s _1 (in the 
2-10 keV range), and may include UCXBs, although the disk 
instability model predicts peak luminosities > 10 37 erg s _1 for 
UCXBs. |King & Wijnands] ( |2006) > found that the luminosities 
of some very faint X-ray transients imply mass transfer rates of 
~ 10" 13 M yr -1 , which is consistent with the behavior of old 
UCXBs. 

The additional angular momentum loss increases the time 
derivative of the orbital period, and as a result the actual number 
of systems at long periods in our prediction based on only gravi- 
tational wave radiation (Figs.[8]and|9]l should be reduced by two 
orders of magnitude. 



4.3. Overprediction of UCXBs with short orbital period 

Since UCXBs with an orbital period shorter than ~ 30 min are 
expected to be persistently bright, our overprediction of these 
binaries (~ 5 - 50 systems based on the disk instability model, 
depending on assumptions in the model, is about one order of 
magnitude more than the three observed Bulge UCXBs) can 
have several causes: the population synthesis model produces 
too many UCXBs, fewer UCXBs survive the onset of mass trans- 
fer, or short-period UCXBs are bright less than 100% of the time. 

It is uncertain whether the white dwarf donor mass limit 
of 0.38 M Q for isotropic re-emission (based on the zero- 
temperature mass-radius relation and a consideration of the en- 
ergy necessary to eject matter) should be used as the threshold 
for the survival of a system. A different assumption in the details 
of the hydrodynamics of the onset of mass transfer could result 
in either a lower or a higher limit. However, the number of sur- 
viving white dwarf donors is not very sensitive to small changes 
in the isotropic re-emission limit because of the small number 
of donors around this value (Fig. [2]). On the other hand, due to 
the sensitivity of this donor mass limit to the actual and criti- 
cal mass transfer rates (van Haafte net al.|201 2b), the number of 
surviving UCXBs could be significantly smaller if the isotropic 
re-emission efficiently is less than unity, which seems plausible. 
However, systems could survive a mass transfer rate exceeding 
the isotropic re-emission limit if an UCXB with a ~ 2 min or- 
bital period is able to survive a few hundred years (~ 10 8 orbits) 
with a significant amount of matter orbiting in and around the bi- 



nary, perhaps in a circumbinary ring ( |Soberman et al.|1997| |Ma| 
& Li |2009 1 that would be removed at a later stage. The value 
of isotropic re-emission limit can be tested only using short- 
period UCXBs, because for those the white dwarf donor channel 
is expected to dominate. In the total population (mostly long- 
period systems), the helium donor channel is more important. 
This channel does not experience the high mass transfer rates 
that characterize the onset of mass transfer from a white dwarf 
donor. Near the period minimum the mass transfer rate remains 
below approximately three times the Eddington limit (Fig. [4) 

Even though Fig.|7]shows that the contribution of the helium- 
star channel dominates, its importance has not been established 
observationally yet, as no detached (short orbital period) helium 
star-neutron star binaries have been discovered so far ( Nelemans 
et al. 2010 1. If helium burning stars would turn into white dwarfs 



before the onset of Roche-lobe overflow, over 90% would be 
unable to survive as a binary system. 

Given these uncertainties, our overprediction by approxi- 
mately one order of magnitude may simply be a consequence 
of poorly known parameters in our simulation, and in this case 
no problematic discrepancy between our results and X-ray ob- 
servations would remain. 



4.4. Consequences for the population of millisecond radio 
pulsars 

If mass transfer would cease completely in most old UCXBs as 
suggested by the existence of PSR J1719— 1438, then a fraction 
of the neutron stars they harbor would become visible to us as 
(binary) millisecond radio pulsars for several billion years. The 
same is probably true for UCXBs in environments closer to us 
than the Bulge. The number of isolated or binary millisecond 
radio pulsars in the Bulge suggested by the number of UCXBs 
with > 70 min periods we predict in our standard model, af- 
ter correcting for the factor of ten overestimation in the UCXB 
birth rate identified above, is about 2 x 10 4 , assuming the pulsars 
have not yet turned off as a result of spinning down. The esti- 
mated Galactic population of mil li second pulsa rs , based on the 
observed population, is 4 x 10 4 ( |Lorimer|2008| . Of the known 
Galactic Disk population (i.e., excluding globular clusters) of 
~ 100 radio pulsars with spin periods shorter than 10 ms, about 
60% have a companion too massive to be consistent with late- 
time UCXB evolution, based on the ATNF Pulsar Catalogu^\m 
January 2013 (Manchest er et al.|20 05). (These have descended 
from hydrogen-rich low-mass and intermediate-mass X-ray bi- 
naries, Bhattacharya & van den Heuvel 1991 Tauris 2011 ) 
Extrapolating this to the estimated Galactic population, ~ 2 x 10 4 
millisecond pulsars are left that have no or a very low-mass com- 
panion, roughly the same number we predict from UCXBs for 
the Bulge alone (after normalizing the number of short-period 
UCXBs). 

It seems likely that millisecond pulsars have evaporated their 
companion entirely and are left as isolated millisecond pulsars 
given the reasonable match between the predicted number of old 
UCXBs and the number of isolated millisecond radio pulsars, 
combined with the very small number of observed millisecond 
radio pulsars with companions with masses lower than 0.01 M . 
Alternative formation channels for isolated millisecond pulsars 
are spin up of a neutron star by a disrupted white dwarf com- 
panion ( |van den Heuvelp984| l, and disruption of a millisecond 
pulsar binary during the supernova explosion of the donor star 
in a high-mass X-ray binary (e.g. Ca milo et al.|19 93). The num- 



http : //www . atnf . csiro . au/people/pulsar/psrcat/ 
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ber of neutron star-white dwarf mergers, however, seems too 
high to be consistent with the number of isolated millisecond 
pulsars that have formed. Based on Fig. [2] merging systems are 
much more common than surviving UCXBs, also after includ- 
ing UCXBs from the helium burning donor channel (see also 
|Iben e t al. 1995). On the other hand, the number of millisecond 
pulsars that lose their companions when it explodes as a super- 
nova seems too small to be responsible for a large fraction of 
the isolated millisecond pulsars (Burgay et al. 2003 1 Belczynski 



eTaT1[20T0 



moreover the high-mass donor star may not live 
long enough to spin up the neutron star to a spin period shorter 
than ~ 10 ms. 



5. Summary and conclusions 

We modeled the present-day population of primordial ultracom- 
pact X-ray binaries in the Galactic Bulge with the purpose of 
gaining insight in their formation and evolution. Both binary 
evolution and accretion physics determine the observable pop- 
ulation, and we attempted to disentangle these in this study. 
We considered three main formation channels: systems that start 
Roche-lobe overflow by a white dwarf donor, a helium burning 
donor or an evolved main sequence donor. Our simulations have 
not produced UCXBs containing a black hole, because most sys- 
tems with a very massive primary merge during unstable mass 
transfer, and the small number that remains is expected to merge 
during the onset of mass transfer from a relatively massive white 
dwarf to the black hole. Thus, all UCXB systems in our simula- 
tions have a neutron star accretor. 

The vast majority of UCXBs form via the helium burning 
donor channel (81%) or the white dwarf donor channel (18%), 
and therefore their exposed cores are expected to show either 
carbon and oxygen in their spectra, or helium, as well as small 
amounts of other reaction products. These two channels differ in 
the delay time between the zero-age main sequence and the onset 
of Roche-lobe overflow to the neutron star. In the white dwarf 
channel this delay can be as long as the age of the Universe or 
more (though for most systems it is less than a few billion years), 
whereas in the helium burning channel the delay is less than 1 
Gyr. 

The size and characteristics of the present-day population 
are only marginally dependent on the assumed width, cr, of the 
Gaussian distribution describing the star formation history if this 
value is < 1 Gyr. This is because these values of cr are small 
compared to the age of a 10 Gyr old system. A broad star forma- 
tion history allows for recent star formation and short orbital pe- 
riod UCXBs with a helium burning donor origin, because of their 
short delay time. With a narrow star formation history, short- 
period UCXBs must have a white dwarf origin and therefore can 
have helium composition. 

Very short period UCXBs can have a helium or carbon- 
oxygen white dwarf donor, since these must have formed re- 
cently. Recent UCXB formation is dominated by the white dwarf 
donor channel, even for cr = 2.5 Gyr. 

The number of predicted systems with orbital periods shorter 
than ~ 30 min is particularly important, since those systems are 
probably observable as persistently bright sources, and therefore 
well suited to test and calibrate the simulations. We predict about 
40 bright sources, mostly of helium and carbon-oxygen compo- 
sition and with orbital periods shorter than 30 min. The UCXBs 
with the shortest periods (< 20 min) are more likely to have 
helium composition. The observed number of bright UCXBs is 
about ten times smaller than suggested by our model, which 
reflects the uncertainties in the adopted star formation history, 



initial binary parameters, natal kick velocities of neutron stars, 
common-envelope parameters and the onset of mass transfer to 
a neutron star accretor. 

We predict about (0.2 - 1.9) x 10 5 UCXBs in the Galactic 
Bulge, and we stress that such a large population is necessary 
based on the simple argument that the evolutionary timescale 
of UCXBs increases rapidly towards longer orbital periods, and 
therefore the observed number of short-period UCXBs, in the 
Bulge and also in the Galactic Disk, implies several orders of 
magnitude more UCXBs at long orbital periods (> 60 min). With 
different model assumptions, this number could be up to an order 
of magnitude lower. 

Irradiation of the donor star by the neutron star and accre- 
tion disk strongly influences UCXB evolution, at least at orbital 
periods longer than 40 min. These systems evolve much faster, 
probably by ~ 100 times, than they would if their evolution was 
driven exclusively by angular momentum loss via gravitational 
wave radiation, as assumed in this paper. UCXBs with orbital 
periods longer than 1 h have not been detected yet, which im- 
plies that, if existent, these systems are very faint in all electro- 
magnetic bands (and therefore cannot be considered true X-ray 
binaries). We suggest that the majority of these systems have or- 
bital periods on the order of 1.5 - 2.5 h rather than the ~ 1.3 h 
expected from gravitational wave driven evolution. Furthermore 
we expect that the neutron stars have companions with masses 
much lower than 0.01 M , and could very well have evaporated 
their companions entirely, being left as isolated millisecond pul- 
sars. 

In a forthcoming paper we will model the population of 
hydrogen-rich low-mass X-ray binaries in the Galactic Bulge. 
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Appendix A: Binary initial mass function and 
normalization of the simulation 

In the initial binary system, the more massive component is 
called the primary. We use primary-constrained pairing to con- 
struct 'zero-age' binaries (Kouwenhoven et al. 2008). The pri- 
mary masses M pl ; mary of the zero-age main sequence binaries are 
drawn from the stellar initial mass function (IMF) of primaries 



in massive star clusters that we derive from the results by Kroupa 
POTH) , where M is the stellar mass and 0.08 < M/M Q < 100. 
The mass ratio < M secon dai- y /M pr i mary < 1 of the components 
is subsequently drawn from a constant distribution ( Krai cheva] 
et al. 1989| Hogeveen T992) - secondary masses lower than 
0.08 Mq are accepted. The eccentricity e distribution is propor- 
tional to e between and 1, and the semi-major axis a distribu- 
tion is invers ely proportional to a ([Popova et al.|1982||Abt|1983[ l, 



up to 10 R Q (Duquennoy & Mayor 1991 1 - the lower limit is set 



by the requirement that the initial stellar radii fit inside the cir- 
cularized orbit. 

The specific binary fraction as a function of M is given by 
the observationally practical definition (Reipur th & Zinnecker| 
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[T993)|Kouwenhoven et al.|2009) 

^binary (-^primary 



S(M) = 



M) 



^single C^O ^binary (^^primary 
AW a ry(M primary = M) 



M) 



(A.l) 



IMF(M) 

where N s i ns \ e (M) is the distribution of single stars of mass M, 
^binary (^primary = M) the distribution of binary systems con- 
taining a primary of mass M, and IMF(M) the IMF of sys- 



tems (single stars and multiple systems combined) by Kroupa 
(|2001| l. Based on observations summarize d in |Kouwenhoven 
|eTaLp|2009) ; |Kraus & Hillenbrand| ( |2009> ; |Sana et al | q20T2] T 
we approximate 



S(M) 
where we 



log 10 (M) (0.08 < M/M < 100) (A.2) 

assume all multiple 
can be separated as 



systems to be binaries. 



(A.3) 



Equation ( A. 1 

AW a ry(M pnmary - M) oc S(M) IMF(M), 

N singie (M) cc [1 - S(M)] IMF(M). 

It follows that single stars are more common than binary sys- 
tems; there are 1.6 single stars for each binary system. The 
mass per binary system including the corresponding single stars 
(which can be a fractional number) is given by 

1 



Mr 



-M 



primary 



"I" ^secondary •> 



(A.4) 



^(^^primary) 

and the average star forming mass for each binary system formed 
(i.e., including mass from single stars) by 

Mj — 1 H I IMF(M)MAM x, 1.9 M Q (A.5) 

Jo.08 Mq \ 



(the factor 1/2 appears because the average secondary mass is 
equal to half of the average primary mass for the chosen con- 
stant mass ratio distribution). This number is the sum of the av- 
erage primary mass (0.86 M B )J^j the average secondary mass 
(0.43 M ) and the corresponding average mass in single stars 
per binary system (0.64 M ). A lower limit of 0.1 M increases 
the average mass per binary by ~ 12%. Overall two-thirds of the 
star-forming mass is in binaries. The total number of binaries 
that forms in the Galactic Bulge is normalized using the total 
number of stars 

-100 M 



5.2 x 10 y 



(A.6) 



f 1UU M ° 1 x 10 1U M G 
N biamy (M)dM = T 2 

Jo.08 M M T 

Of all primaries, 1.3% have a mass higher than 8 M . For these 
masses, the power-law slope of the primary IMF (defined over 
linear mass intervals), from which we draw primary masses, 
varies between -2.15 (for M = 8 M Q ) and -2.2 (M = 100 M Q ), 



compared to the estimate of -2.3 by Kroupa (2001 1 for the com- 
bined IMF of single stars and primary components. The IMF of 
primary components Aginary (Mprim ai-y = M) is flatter than the IMF 
of systems IMF(M) because Eq. (A.2i is an increasing function 



(most low-mass stars are single whereas massive stars are usu- 
ally in binaries)[^]The IMF for single stars only is steeper than 
-2.3 and steepens towards high mass. 

11 This value is higher than the average mass of 0.57 M Q of the IMF 
by Kroupa (20011 because single stars are excluded. 



12 A non-zero mass-independent binary fraction leads to an IMF of all 
stars combined that is steeper than the IMF of systems ( Sagar & Richtler 



[T99T] Scalo 1998 , Kroupa 2001 ). This does not affect our method as we 
only consider the IMF of all primary components. 
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